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The properties of two measures of charge fluctuations D and A$ 9 are discussed within 
several toy models of nuclear collisions. In particular their dependence on mean particle 
multiplicity, multiplicity fluctuations and net electric charge are studied. It is shown 
that the measure A$ 9 is less sensitive to these trivial biasing effects than the originally 
proposed measure D. Furthermore the influence of resonance decay kinematics is analysed 
and it is shown that it is likely to shadow a possible reduction of fluctuations due to QGP 
creation. 
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I. INTRODUCTION 



Recently Jeon, Koch M and Asakawa, Heinz and Mueller Q| argued that the study of event-by-event 
fluctuations of electric charge in high energy nucleus- nucleus collisions may provide information on the 
state of matter in an early stage of the collision. The authors calculated the magnitude of the charge 
fluctuations in the Quark-Gluon Plasma (QGP) and in a hadron gas. These calculations show that the 
fluctuations in the QGP should be significantly smaller than in a hadron gas. Thus they can be used as a 
signal for the creation of the deconfined phase providing that the initial fluctuations survive hadronization 
and their relaxation time is significantly longer than the time of the hadronic stage of the collisions |l| || . 

The charge fluctuations should also be sensitive to other effects like the number of resonances at 
chemical freeze-out [|]j5| and fluctuations occurring in the initial stage |6) . Consequently their analysis is 
interesting also beyond the QGP hypothesis. 

First experimental results on charge fluctuations in central collisions of heavy nucleus at SPS and 
RHIC || are already available and don't show the expected reduction of fluctuations due to QGP creation. 

The analysis of charge fluctuations is relatively easy to perform experimentally due to the typically 
good resolution of the measurement of electric charge in tracking detectors positioned in a magnetic field. 

The magnitude of the measured charge fluctuations is, however, also dependent on trivial, uninteresting 
effects, which may shadow the studied physics. The two most important of these effects are: 

- the fluctuations in the event multiplicity, introduced mostly by the variation of the impact parameter; 

- changes in the mean multiplicity due to changes of the colliding system (e.g. p+p, Pb+Pb), changes 
of collision energy and due to changes of the size of the acceptance in which fluctuations are studied. 

In order to minimize the sensitivity to these effects two measures of fluctuations were proposed: D H 



and $ p 11 



The aim of this paper is to study the properties of these measures using simple toy models of particle 
production in nuclear collisions. 

In Section II the fluctuation measures will be introduced. Their properties will be tested in Section 
III. Summary and conclusions are given in Section IV. 

II. MEASURES OF CHARGE FLUCTUATIONS: D AND A$ Q 

The measure D is defined as follows pl^| : 
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N pos and N neg are measured multiplicities of positively and negatively charged particles within the 
analysed acceptance. The symbol (...) represents averaging over events. (Nch) is the mean number of 
charged particles within the acceptance and (N c h) tot denotes the mean total event multiplicity of charged 
particles in full phase space. The factors C y and were introduced to remove the influence of global 
charge conservation (C y ) and effect of non-zero net charge (C^). Under the assumption of Poissonian 
distribution of N pos and N neg (providing the number of accepted particles is much smaller than the total 
number of particles =4> C y ~ 1) and assuming zero net charge (Q = N poStto t — N negito t = = 1) one 
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obtains D = 4 |1]J9|. Many hadrons measured in the final state originate from the decay of resonances. 
Therefore those are correlated in kinematic quantities like rapidity. Consequently if we consider an 
acceptance window in a kinematic variable which is large compared to the mean separation of the decay 
products in this variable the probability that all decay products will fall into it is high and therefore 
the charge fluctuations should be reduced significantly ]j],||J|- I n the QGP phase the unit of charge is | 
instead of 1 as in the hadron phase. Smaller charge units cause smaller charge fluctuations. Assuming 
that the QGP fluctuations are frozen one expects JjMl' 

A well established measure of event-by-event fluctuations is the variable $ [[l0| which is defined in the 
following way: 



where: 



N 



z = x — x, Z = ^^(xi — x). (7) 

a; is a single particle variable, N is the event multiplicity within acceptance and overline denotes averaging 
over a single particle inclusive distribution. By construction, for a system which is an independent sum of 
identical particle sources the value of <& is equal to the value of $ for a single particle source independent of 
the number of superimposed particle sources and its distribution . In the original proposal JlO| <E> was 
introduced as a measure of fluctuations of a kinematic variable like transverse momentum. Furthermore 
it was proposed jflj to use $ as a measure of fluctuations of discrete quantities. Following this suggestion 
in this publication x is identified with the particle charge q. Some properties of the $ measure for discrete 
quantities were derived in Ref. . For a scenario in which particles are correlated only by global charge 
conservation (gec) and providing Q — we calculated the value of $ g : 



$ q ,gcc = (8) 

where 

P=7pf- (9) 

{-l^ch/tot 

Thus the value of & q ,g C c depends only on the fraction of accepted particles. In order to remove the 
influence of gec we suggest therefore to use a difference: 

A$ g = - $ 9 , scc . (10) 

Consequently the value of A$ 9 is zero if the particles are correlated only by global charge conservation. 
It is negative in case of an additional correlation between positively and negatively charged particles and 
positive if there are anticorrelations. One possible source of additional correlations could be the creation of 
a QGP. The first part in the definition of $ takes into account all, statistical and dynamical fluctuations. 
The second part quantifies only the statistical ones. Therefore the value of $ calculated for different 
systems is independent of their degrees of freedom (quarks and gluons or hadrons) of the considered 
system if its fluctuations are only statistical. However, experimentally we measure the fluctuations 
of particles only after hadronization. A hadronization which conserves charge and entropy produces 
correlations. Consequently the measured value of $ g and A<j> g should be smaller, if a QGP was created 
in the early stage of the collision. 



III. TOY MODELS OF NUCLEAR COLLISIONS 



In this section properties of the fluctuation measures, D and A$ 9 , are discussed. First we will assume 
that the only source of correlations is the global charge conservation. Under this assumption we will study 
the behaviour of the measures at small (N^) and their response to non-zero net charge and fluctuations 



of event multiplicity. Furthermore we also calculate their values in two toy models, a 'p'-gas model and a 
QGP model in which additional correlations are present. In the study we use a simple Monte Carlo (MC) 
procedure, which simulates the particle distribution in momentum space. The values of D and A$ 9 are 
calculated for particles falling into rapidity windows centred around mid-rapidity. In the MC simulations 
we used values of total multiplicity ((-/V c /j) tot ) and net charge (Q) which are typical for central Pb + 
Pb collisions at 40 A-GeV as measured by the NA49 collaboration Q]. In the figures the measures are 
plotted as a function of (N c h), which is a positive monotonic function of the size of the rapidity window. 
The ranges of the measures D and A<5> 9 plotted on the vertical scales are chosen in order to make both 
measures similar sensitive to resonance decay kinematics and QGP creation. Statistical errors of the MC 
simulations are plotted but in most of the cases they are smaller than the symbol size. 

A. Model A 

In the first model we analyse the dependence of the measures D and A<i> 9 on (N c h). It is assumed 
that the total multiplicity is fixed, N c h^tot — 850, and the net charge in each event is zero (Q = 0). The 
dependence of D and A$ g on (N c h) is shown in Fig. 1. 




FIG. 1. The dependence of D and A$ q on (N c h) obtained in model A. 



As expected the value of A$ g is zero and the value of D is close to 4 for large (N c h). However for (N c h) 
smaller than 200 a significant deviation of D from 4 is observed. It is because the assumptions used to 
derive the value D — 4 are not fulfilled in the limit (N c h) — ► 0. 

B. Model B 

Both measures are constructed in order to be independent of the fluctuations of N^tot- To test it we 
varied the total number of charged particles from event to event according to a Gaussian distribution 
with a mean value 850 and a dispersion of 60 and 120. As in the model A, the global net charge was 
assumed to be zero. The results of this simulation are shown in Fig. 2. 
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FIG. 2. The dependence of D and A$ 9 on (N c h) and a{N c h,tot) obtained in model B. 

The behaviour of D and A$ g is similar to this obtained in Model A. Thus the influence of multiplicity 
fluctuations is small, especially on A$ g . 

C. Model C 

After we have observed that the measures are insensitive to fluctuations of N c h,tot we turn to check 
the influence of non-zero net charge. We assume N c h } tot = 850 and Q — 60. Fig. 3 shows the results. 
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FIG. 3. The dependence of D and A«I> 9 on (N c h) and net charge Q obtained in model C. 

The measure D is explicitly corrected for non-zero net charge by the factor C M , but it still shows a 
significant dependence on it. Much weaker dependence on the net charge is observed for A$ g . For an 
acceptance larger than 50% its value is somewhat higher than the calculated one assuming Q = 0. 

D. Model D 

In nuclear collisions the number of charged particles fluctuates and the net charge is approximately 
proportional to N c h,tot- To see how the measures behave under these assumptions we used a model 



in which N c h,tot was generated from a Gaussian distribution {{N c h,tot) — 850, <y{N c h,tot) = 80) and 
Np OS ,tot = j2 ■ N c h,tot, N neg ,tot = j2 ■ N ch,tot • The results of this simulation are presented in Fig. 4. 




FIG. 4. The dependence of D and A$ 5 on (N c h), for fluctuating N c h,tot and for net charge which is proportional 
to N c h,tot obtained in model D. 

The value of D is significantly higher than 4 in the full range of {N c h) ■ A$ 9 is much less dependent on 
the studied effects. It deviates slightly from for (N c h) > 0.5 ■ (iV c /, fo$). 



E. QGP Model 

Following Ref. jjj we construct a simple toy model of the QGP. Under the assumption of zero baryo- 
chemical potential and providing zero quark mass we calculate the ratio of up-, down-, antiup- and 
antidown- quarks and gluons in equilibrium. Assuming entropy and net charge conservation in each 
rapidity window during the evolution from the QGP to the final hadron state the number of pions (N) 
and their net charge is calculated. The number of charged pions is taken to be N c h = | • N based on 
isospin symmetry. Using this model we obtain the results shown in Fig. 5. 




FIG. 5. The dependence of D and A$ 9 on {N c h) obtained in a QGP toy model. 
Our simulation confirms results published in Ref. B, the D value in our QGP toy model is close to 1. 



The new result is the calculation of A<I> g in this model. Its value is mostly significantly smaller than zero. 
It increases from —0.5 for small (N c h) to —0.1 for the largest multiplicity. In case all particles are used, 
there is no net charge fluctuations due to charge conservation and therefore the value of A$ g should be 
0. We note that the previously discussed trivial effects were increasing D and A$ g to values above 4 and 
0, respectively. Thus they can not mimic suppression of fluctuations observed in the QGP toy model. 



F. 'p'-gas Model 

In this model we study the influence of the resonance decay kinematics on D and A$ g . We use a MC 
simulation which generates 400 neutral resonances ('p') in each event. The mass of the 'p' is fixed to 
770 MeV/c 2 . The 'p' rapidity distribution is Gaussian with a(y) = 0.8, 1, 2 or 3 and the transverse mass 
distribution is 'thermal': dr f" d = C ■ m T ■ e~ mT/T with T = 170 MeV. All 'p's decay into two charged 

pions which are then used to calculate values of D and A<i> g . In Fig. 6 we show the results of a simulation 
for cr(y) = 1.0. 
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FIG. 6. Dependence of of D and A3> 9 on {N c h) obtained in 'p'-gas model for o(y) = 1.0. 



For low (N c h) the values of D and A<& 9 are higher then 4 and 0, respectively. With increasing (N c h) 
these values decrease. Resonance decay kinematics influences charge fluctuations in two different ways. It 
dilutes the effect of global charge conservation if only one of the decay products falls into the acceptance. 
This will increase the values of D and A& q . If both decay products fall into the acceptance (N c h) increases 
but the net charge doesn't change. D and A<I> 9 measure charge fluctuations normalised by the number 
of charged particles. Therefore their values will drop in this case. Which effect dominates depends on 
the size of acceptance window (Ay, upper scale on the plots) , mass of the decaying resonances and their 
rapidity distribution. The sensitivity on the last one is shown in Fig. 7. 
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FIG. 7. Dependence of D and A$ 9 on (N c h) obtained in 'p'-gas model for cr(y) = 0.8, 2.0 and 3. Note that 
horizontal upper scale (Ay) changes with ff(y). 

At SPS energies a(y) of the p distribution is expected to be close to 1 In Figs. 6 and 7 (upper 

plots) we see that small changes of <r(y) around 1 have strong influence on the values of D and A<I> g . 
For er(y) = 1 we observe a transition from fluctuation enhancement to suppression with an increasing 
rapidity window. For <r(y) = 0.8 the values of D and A<j> 9 are above or equal to their values expected for 



independent particle production with global charge conservation for all Ay. We would like to stress here 
that in our model the 'p's are neutral and therefore there are no charge fluctuations before the 'p's decay. 
In Fig. 6 and Fig. 7 (upper plots) we see that for cr(y) values expected at CERN SPS the resonance decay 
kinematics increases the initial zero charge fluctuations to values which are close to 4 or 0, respectively. 
The fluctuations are significantly lower for larger values of er(y), e.g. er(y) = 2 and 3, as are expected at 
RHIC and LHC energies. In the case of QGP creation the measured fluctuations will be composed of the 
initial QGP fluctuations and fluctuations produced by the decay of resonances. The latter ones are, as 
our model shows, still larger than the predicted QGP fluctuations even for cr(y) = 2 and 3. Therefore we 
conclude that resonance decay kinematics may significantly shadow charge fluctuations developed at the 
early stage of the collision. 



IV. SUMMARY AND CONCLUSIONS 



The properties of two measures of charge fluctuations D and A<j> g were studied within several models. 
We have shown the dependence of both measures on trivial effects, in particular, on number of particles 
used per event, on the net electric charge and on fluctuations in the total event multiplicity. Furthermore 
their response to resonances decay kinematics and to the reduction of fluctuations due to QGP creation 
were analysed, too. From our simulations we learn that the measure A$ g is significantly less dependent 
on trivial effects than D and that both measures are sensitive to frozen QGP fluctuations and to resonance 
decay kinematics. Therefore we suggest to use A$ g as a proper measure of charge fluctuations. Using 
the 'p' gas model we have shown that the influence of resonance decay strongly depends on their rapidity 
distribution and on the acceptance window. Given the typical widths of the rapidity distributions of p's 
at CERN SPS energies they can easily increase charge fluctuations from QGP values to values expected 
for independent particle production with global charge conservation. A further development of methods 
which allow to distinguish various sources of charge fluctuations is obviously necessary. 
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